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Abstract 



In this paper, we consider matrices whose inverses are tridiagonal Z-matrices. Based on a 
characterization of symmetric tridiagonal matrices by Gantmacher and Krein, we show that a 
matrix is the inverse of a tridiagonal Z-matrix if and only if, up to a positive scaling of the rows, 
it is the Hadamard product of a so called weak type D matrix and a flipped weak type D matrix 
whose parameters satisfy certain quadratic conditions. We predict from these parameters to 
which class of Z-matrices the inverse belongs to. In particular, we give a characterization of 
inverse tridiagonal M-matrices. Moreover, we characterize inverses of tridiagonal M-matrices 
that satisfy certain row sum criteria. This leads to the cyclopses that are matrices constructed 
from type D and flipped type D matrices. We establish some properties of the cyclopses and 
provide explicit formulae for the entries of the inverse of a nonsingular cyclops. We also show 
that the cyclopses are the only generalized ultrametric matrices whose inverses are tridiagonal. 
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1 Introduction 



In many mathematical problems, Z-matrices and M-matrices play an important role. It is 
often useful to know the properties of their inverses, in particular when the Z-matrices and the 
M-matrices have a special combinatorial structure. In this paper, we investigate the properties 
of inverse tridiagonal Z-matrices and M-matrices, i.e., matrices whose inverses are tridiagonal 
Z-matrices or M-matrices. We also highlight some connections between weak type D matrices 
(a class that generalizes type D matrices as defined by Markham and inverse tridiagonal 
Z-matrices. 

First, under the assumption of irreducibility, we show that a matrix is the inverse of a tridiag- 
onal Z-matrix if and only if, up to a positive scaling of the rows, it is the Hadamard product 
of a weak type D matrix and a flipped weak type D matrix whose parameters satisfy certain 
quadratic conditions (Theorem . This characterization parallels (and is based on) the char- 
acterization of (symmetric) Green matrices by Gantmacher and Krein Further, recalling 
the classification of Z-matrices by Fiedler and Markham |^] , we predict the class Lg of a tridi- 
agonal Z-matrix based on the parameters of the associated weak type D matrices (Theorem 



3.4). In particular, we find conditions on the parameters so that the inverse is a tridiagonal 



M-matrix (Corollary p. 61) . 

Next, we associate type D matrices with tridiagonal Z-matrices via the so called cyclopses. 
These are matrices that admit a block partition comprising two diagonal blocks that are of 
flipped type D and of type D, respectively, and two off-diagonal blocks that have constant 
entries. We find conditions on the parameters of the associated type D matrices and the 
constant off-diagonal entries so that the inverse of a cyclops exists and is a tridiagonal Z- 
matrix; its nonzero entries are also found explicitly in terms of the parameters (Theorem [4.6| ). 
When a cyclops is a priori nonsingular, we provide necessary and sufficient conditions so that 



its inverse is a tridiagonal Z-matrix (Corollary 4.7); as before we can predict the class Lg of 
the tridiagonal Z-matrix (Theorem |4.8|). 



Cyclopses (with nonnegative entries) were encountered by the authors as a special case of the 
generalized ultrametric matrices (see [^] and [^]), which is a class of inverse (row and column 
diagonally dominant) M-matrices. We conclude by finding necessary and sufficient conditions 
so that a cyclops is the inverse of a (row and column) diagonally dominant tridiagonal M- 
matrix (or equivalently a totally nonnegative generalized ultrametric matrix) (see Theorems 
4.10 , [4.12| , and 4.14| ). These results amount to a characterization of the generalized ultrametric 



matrices whose inverses are tridiagonal. 

We continue with the precise definitions of the terms mentioned above and the notational 
conventions. 
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2 Preliminaries 



We let e denote the all ones vector and ej the j-th standard basis vector in ]R^. Given a positive 
integer n we let (n) = {1,2, .. . ,n}. Let o denote the Hadamard (i.e., entrywise) product of 
matrices. For A = [uij] E IR°''^, by A{i\j) we denote the submatrix of A obtained by deleting 
the i-th row and the j'-th column. Given R,SC (n) we write Ajis for the submatrix of A 
whose rows and columns are indexed by R and S, respectively. If 5 = (n) \ i? and if Ajiji is 
nonsingular, then the Schur complement of ^_r_r in A is defined and denoted by 

A/Arr = Ass - Asr{Arr)-^Ars- 
It is well known that det^ = detARRdet{A/ARR). 

We call A = [uij] £ IR'^'" a Z-matrix if mj < for all i ^ j. For any nonnegative integer 
s < n we denote by Ls the set of all matrices A = tl — B £ IR'^'", where B is an entrywise 
nonnegative matrix and where Ps{B) < t < Here Ps{B) denotes the maximum among 

the spectral radii of all the s x s principal submatrices of B (we take po = —oo and Pn+i = oo)- 
In particular, A is an M~matrix if it can be written as ^ = tl — -B, where B is an entrywise 
nonnegative matrix and p{B) := Pn{B) < t. 

The next theorem, found in and is a characterization of the nonsingular Z-matrices 
in Lg. 

Theorem 2.1 Let A € M'^'" be nonsingular Z-matrix. Then A Lg if and only if one of the 
following alternative cases a) or b) holds: 

a) (i) detA < 0, 

(ii) all principal minors of A^^ of order greater than or equal to n — s are nonpositive, 
and 

(iii) there exists a positive principal minor of A~^ of order n — s — 1. 

b) (i) det^ > 0, 

(ii) all principal minors of A~^ of order greater than or equal to n — s are nonnegative, 
and 

(iii) there exists a negative principal minor of A~^ of order n — s — 1. 

Markham defined in |^ type D matrices as follows: A = [aij] G IR"'"^ is of type D (or a type D 
matrix) if 

{Oi, i < j, , 
where an > a„_i > . . . > ai. 
aj, i> j, 

We refer to the ai (i = 1, 2, . . . , n) as the parameters of A. We also consider similarly con- 
structed matrices, without constraints on the parameters Oj, to which we refer as of weak type 
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D. Moreover, we call A a flipped type D matrix (resp., a flipped weak type D matrix) if PAP^ 
is a type D matrix (resp., a weak type D matrix), where P is the permutation that reverses the 
order of the indices 1,2, ... ,n. We enumerate the parameters of a weak type D matrix, as well 
as the parameters of a flipped weak type D matrix in a way such that the i-th parameter is 
equal to the i-th diagonal entry of the matrix. To illustrate these definitions and the relevant 
notation, let 





"-1 


-1 


-1" 




"-3 


2 


1" 


A = 


-1 


2 


2 


, B = 


2 


2 


1 




-1 


2 


3 




1 


1 


1 



Then A is of type D with parameters Oj given by (—1,2,3) and B is of flipped weak type D 
with parameters bi given by (—3, 2, 1). 

Gantmacher and Krein defined in Q a Green matrix to be a matrix G G H"^'" such that 
G = A o B, where ^ is a weak type D matrix, i? is a flipped weak type D matrix. The name 
Green matrix is not the only name for these matrices. Originally Gantmacher and Krein called 
such matrices einpaarig or matrix of a couple. Moreover, Markham defined the type D matrices 
as a special case of Green matrices. 

In our discussion, we shall also refer to the following matrices that were introduced in |jl3[ and 
||9|. We say C = [cij] € IR'^''^ is a generalized ultrametric matrix if 

(i) C is entrywise nonnegative, 

(ii) Cii > max{cy,Cji} for ah i,j G (n), 

(iii) every subset of (n) with three distinct elements has a labeling {i,j. A;} such that 

(a.) Cij — Cjfc, 

(b) Cji — Cfcj, 

(c) mm.{cjk,Ckj} > mm{cji,Cij}, 

(d) max{cjfc,Cfej} > maxjcji, Qj}. 



In the aforementioned papers, it is shown that if a generalized ultrametric matrix is nonsingular 
then its inverse is a row and column diagonally dominant M-matrix. 

Next, we introduce a class of matrices constructed from type D matrices; as we show in Section 
4, it contains matrices that are under certain additional conditions are inverse tridiagonal Z- 
matrices. Let C S IR'^''^ and let m < n be a nonnegative integer. We call C a cyclops with eye 
m+ if 



(2.1) 



C 



Cii biEi2 

62 -£-21 C22 
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where Cn is a m x m flipped type D matrix and C22 is a (n — m) x (n — m) type D matrix, 
viz., 



C 



11 



ai 
02 



02 
02 



ar 



ar. 



O'm—1 Q-m 
O^m— 1 Q-m 



22 



Om+1 am+2 



Om+1 Om+2 
Om+1 Om+2 



Clm+2 Om+2 



'^n— 1 Oifi 



with 
(2.2) 



Oi > 02 > . . . > flm and Qn > Qn-l > • • • > O. 



and where i?i2 and E21 are ah ones matrices of appropriate sizes. We refer to the Oj (i = 
1, 2, . . . , n) and 61, 62 as the parameters of the cyclops C. 

In the remainder of this paper, when we refer to a type D matrix, a weak type D matrix, a 
Green matrix, or a cyclops, we assume that the reader recalls the notation and the associated 
parameters indicated in this section. 



3 Hadamard Products of weak type D Matrices 



Gantmacher and Krein proved the following results. 



Theorem 3.1 (Gantmacher and Krein Q) Let G G IR'^''^ he symmetric. Then the following 
are equivalent: 

(i) G is a nonsingular Green matrix. 

(ii) G^^ is an irreducible tridiagonal matrix. 



Lemma 3.2 (Gantmacher and Krein Q) Let G S JR^'"^ be a Green matrix with associated 
parameters Oj, bi. Let hi := aibi-i — ai^ibi for i = 2, 3, . . . , n. Then 

n 

detG = aibnY\_hi. 

i=2 

Moreover, 



detG{i\j) 



detG/hi+i if\i-j\ = l 
if \i-j\ > 1. 
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We proceed by characterizing inverse tridiagonal Z-matrices in the spirit of Theorem 



3.1. 



Theorem 3.3 Let C € M"^'" be nonsingular and irreducible. Then the following are equivalent: 



(i) C ^ is a tridiagonal Z-matrix. 

(ii) There exists a positive diagonal matrix D 

weak type D with parameters ai, and B is 
that aibi-i — Oi-ibi > for all i = 2,3, . . 



€ W'^ such that DC = Ao B, where A is of 
of flipped weak type D with parameters bi, such 
.,n. 



Proof: Let hi := Uibi-i — Oj-i^j for i = 2,3, . . . ,n and suppose that (i) holds. As is well known 
(see 1^ and flj]), there exists a positive diagonal matrix such that C^^D~^ is symmetric. 



Thus, by Theorem 3A, DC is a Green matrix. We also have that C ^D ^ = [ya], where 



-iy+^det{{DC){j\i))/det{DC). 



Hence, by Lemma |3.2| , the superdiagonal entries are 

-1 



(3.3) 



h 



(i = l,2,...,n). 



Thus all hi are positive. Conversely, if (ii) holds, ( ^.3| ) and Theorem 3J imply that C ^D ^ 
is a tridiagonal Z-matrix. Hence is a tridiagonal Z-matrix. ^ 

In the following theorem we determine the class Lg to which an inverse tridiagonal Z-matrix 
belongs. 



Theorem 3.4 Let C € ]R°'" be an irreducible inverse tridiagonal Z-matrix. Let D, A, and B 
be as in condition (ii) of Theorem \3. 4 Then the following hold: 

(i) // detC < then C~^ G Lg with s = mm{t -2,n-r - 2}; 

(ii) if detC > then C^^ G Lg with s = min{t — 2,n — q — 2}, 



where 

^ _ \ n + 2 if aibi+j — ai+jbi > for all i,j 

1 min{j G (n)| there exists an i ^ in) with Oibi^j — aij^jbi < 0} otherwise, 

^ _ f -1 if Oibi+j <0 for all i,j 

1 max{j G {n)\ there exists an i £ (n) with Oibi^j > 0} otherwise, 

_ f -2, if aibi+j > for all i,j 

1 max{j G there exists an i £ (n) with aibij^j < 0}, otherwise. 
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Proof: For i,j € (n) with i > j, define hij = aihj — ajbi. Since C is an inverse tridiagonal 
Z— matrix we liave, by Theorem |3.3|, that > for all i = 2, 3, . . . , n. Moreover, by Lemma 



3.2, we have that 



(3.4) det(Z)C) =ai6„n^M-i> 

i=2 

where D is the positive diagonal matrix chosen in the proof of Theorem We will proceed 
by considering the signs of the principal minors of DC and by applying Theorem Since 
principal submatrices of Green matrices are also Green matrices, the principal minors of DC 
are given by formulae similar to (|3.4| ), and their signs are determined by the corresponding 
quantities hij and and bi. 

First, suppose detC < 0, i.e., ai6„ < 0. Without loss of generality, we can assume that oi > 
and bn < 0. When t ^ n + 2, the definition of t and (3.4) imply that there exists a principal 



submatrix of order n — (t — 1) of DC with positive determinant. (This principal submatrix is 
obtained by deleting rows and columns z + + j — 1, where i,j are the minimal indices 

in the definition of t.) For all principal submatrices of order greater than n — (t — 1), the 



relevant /ijj appearing in the determinantal formula of Lemma 3.2 are positive. It is also clear 
from the definition of r that there exists an (r + 1) x (r + 1) principal submatrix with positive 
determinant. Moreover, all principal submatrices of order h with n > r + 1 satisfy aihn < 0, 
where di and ba are the corresponding parameters. ^From these cases, we obtain the following: 
if s = min{t — 2, n — r — 2}, then there exists a principal minor of order n — s — 1 that is positive. 
Also, all principal minors of order greater than n — s — 1 are nonpositive. Thus, by Theorem 
2.1, C^^ € Lg, showing (i). Similarly we obtain (ii). g 



It is shown in |11] that if C ^ £ Lg and [^J < s < n, then detC < . For inverse tridiagonal Z- 
matrices, this result can be established by considering the changes of the signs in the sequences 



of the parameter Oj and 6j. With q as in Theorem 3.4, if detC > and C is not entrywise 
nonnegative, one obtains that q + 1 > [^J . 



Example 3.5 In the following examples we apply Theorem 3.4 
(i) Consider the matrix 



C 



-2 -2 

-1 -1 

2 2 

2 3 



for which we can write DC = A o B, where D = I, A is of type D with parameters a, given 
by (-2,-1,2,3), and B is of (flipped) weak type D with parameters bi given by (1,1,1,1). 
Notice that Uibi-i — ai-ibi > for all i = 2,3, 4. Moreover, detC < 0, t = 6 and r = 1. Thus 
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(ii) Let 



C = 



-24 


-20 


-16 


-4 




-4 


-4 


-4 


-4 




6 


5 


4 


1 


-20 


-10 


-8 


-2 




-4 


-2 


-2 


-2 




5 


5 


4 


1 


-16 


-8 


16 


4 




-4 


-2 


4 


4 


o 


4 


4 


4 


1 


-4 


-2 


4 


5 




-4 


-2 


4 


5 




1 


1 


1 


1 



C is the Hadamard product of a weak type D and a flipped weak type D matrix with parameters 
(—4, —2, 4, 5) and (6, 5, 4, 1), respectively. By Theorem |3.3| , C^^ is a tridiagonal Z-matrix and 
detC < 0. Moreover, r = 1 and t = 2 and hence S Lq. 



The results above yield the following characterization of inverse tridiagonal M-matrices. 



Corollary 3.6 Let C € 'M^'^ he nonsingular and irreducible. Then the following are equiva- 
lent: 



(i) C is a tridiagonal M-matrix 

(ii) There exists a positive diagonal matrix D G H'^'" such that DC = A o B, where A is of 

a weak type D with parameters Oj, and B is of flipped weak type D with parameters hi, 
such that all the parameters have the same sign and 

0<^<^<...<^. 

Ol 02 On 



Proof: 



(i) implies (ii): As C G L„, we have that detC > and that C is a tridiagonal Z-matrix, 
The implication now follows from Theorems 3.2 and 3.4. 



(ii) implies (i): By Theorem 3.3, C ^ is a tridiagonal Z-matrix. Since C is entrywise positive 
and since every inverse positive Z-matrix is an M-matrix, (i) holds. ^ 



4 Cyclopses 



In this section, we consider inverse tridiagonal Z-matrices that satisfy certain row sum and 
column sum criteria. This leads to a new class of matrices that we have defined as cyclopses 
in Section 2. We begin with some auxiliary results. 



Lemma 4.1 ([0, Observation 3.8]) Let A G M'^'" he of type D with parameters ai and such 
that ai 7^ 0. Then 

A-^e= — ,0,0,...,0 
ai 
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^From Lemma 3.2, for a type D matrix with parameters a^, we obtain that 

n 

(4.5) detC = ai - flj-i). 

j=2 



Lemma 4.2 Let C G H'^'" 6e a cyclops with eye m+ parameters Oj, 6i, 62- Suppose that 7^ 
and a^rt+i 7^ 0. Then 

rir -r ^^^^ jr rir -r 



where Ei, E2 are all ones matrices of appropriate sizes. 



Proof: Follows from Lemma iA and the fact that C /Cu = C22 — &ife2-E'2iC'i/^: 



12- 



Theorem 4.3 Let C G IR^'^ be a cyclops with eye m+ and parameters 04,61,62. Then 

m— 1 n 

detC = (amOm+i - 6162) J| {am-j - am-j+i) {aj - aj-i). 

j=l j=m+2 



Proof: 



Case I (a-m 7^ or am+i 7^ 0): if am 7^ 0, then by Lemma 4.5, Cu is nonsingular and 

detC = detCiidet(C7/Cii) 



m—l 



) det(C22 - 6162/am^i) 



j=l 



m—l 



Oim ifi"m—j 



)(am+l 



6162 



m—l 



6162 



)) 



j=m+2 



= (Omflm+l - 6162) Jl (flm-j " Om-i+l) H ~ "^J-l)- 

j=l j=m+2 

If 

flm+1 7^ the result follows in a similar manner. 

Case II (om = flm+i =0): if = cim+i = and either 61 = or 62 = 0, then C has a row of 
zeros (and thus zero determinant) and the result follows. Assume that 61 7^ and 62 7^ 0. 
Let i? = {m, m + 1}, 5 = {1, 2, . . . , m - 1}, T = {m + 2, m + 3, . . . , n} and [/ = S" U T. 
Then Crr is nonsingular since Om-i > = 0. Hence 



detC = detCRi?det(C/CRj?) = -6162 det 



Css 
Ctt 
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m— 1 n 

— 6l62(am-iam+2) JJ^ («m-j — Om-i+l) Yl i^j ~ ^j-l) 
j=2 j=m+3 

m—1 n 
j=l j=m+2 



The following is an immediate consequence of the above theorem. 

Corollary 4.4 Let C E ]R'^'° he a cyclops with eye m+ and parameters a,, 61, 62. Then C is 
nonsingular if and only if amam+i — ^1^2 7^ 0. Moreover, sgn(detC) = sgn(amam+i — ^1^2)- 

Next, we shah give explicit formulae for the entries of the inverse of a nonsingular cyclops. We 
first need another result on type D matrices proved in [Q. We denote by <8> the Kronecker 
product of matrices. 

Theorem 4.5 Let A be a nonsingular matrix of type D with parameters Oj. Then the inverse 
of A is given by 

n 

(4.6) A-' = (a, - a,_i)-\ 

1=1 



(4.7) vf = { I if j = i 



with oq = 0. Here the vectors u^*) = [vf'] are defined as 

-l if j = i - I 
I if j = i 
otherwise. 

It follows that the entries aij of the inverse of a type D matrix A G IR'^''^ are zero except for 

1 1 1 

an = \ , Onn ~ 



ai 02 — ai an — ffln-i 

11/ 

aii = \ (i = 2,3, ... ,n- 1), 

Oj — Oj-i ttj+i — ai 

1 



O.i+1 — Oi 

Similarly, for a flipped type D matrix A we have 

1 1 1 

an = , Onn = \ , 

Oi — 02 a„ fln-l — CLn 

aii = \ (i = 2,3, ... ,n- 1), 

o-i — aj+i aj_i — ai 

1 

oti^i+i — ai+i,i 



o,i — aj+i 
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Theorem 4.6 Let C E ]R°''^ be a cy clops with eye m+ and parameters a^, 61, &2- Suppose that 
QmO'm+i — ^1^2 7^ 0. Then A = = [aij] exists and is a tridiagonal matrix with entries 
given by 

1 

an = , 

ai — 02 



Moreover, 



an 



^ + (z = 2,3,...,m-l), 



ar, 

0^771+1,771+1 

1 



+ 



1 



1 



'^TTl 



(i = l,2,...,m-l) 



+ 



1 



an 



+ 



0"m,am+l — blh2 am+2 — CLm+1 
1 



«i — Oj-l Oj+l — Oi 



(i = m + 2, m + 3, . . . , n — 1) 



1 



Or- 



1 



«i,i+l — Oj+l,' 



Om m+1 



bi 



(i = m + l,m + 2,...,n — 1). 

b2 



u u ' Om+1,771 — u u ' 

amO-m+l — O1O2 amam+1 — O1O2 



Proof: Recall that our assumption that amCim+i — ^1^2 7^ is equivalent to C being invertible. 



Case I (am 7^ and a^+i 7^ 0): if Om 7^ then Cn is a nonsingular flipped type D matrix 



and hence it is the inverse of a tridiagonal Z-matrix (by the results in [11 1). If Om+i 7^ 0, 
then C22 is a nonsingular type D matrix and hence it is also the inverse of a tridiagonal Z- 



matrix. By Lemma C/Cu is of type D and C/C22 is of flipped type D. Moreover, by 
( [4.5[ ) applied to C/Cn and C/C22 and since amflm+i — 6162 7^ 0, both Schur complements 
are nonsingular and thus (using formulas from P, (10), p. 773]) 



A 



(0/022)-' 

-b2022 E2l{0 / O22)'' 



-biC^^'Eu(C/C-_ 



11; 



^-1 



(C/C: 



11; 



Since O/O22 is of flipped type D, its inverse is a tridiagonal Z-matrix with all row sums 
zero except the m-th (last). Since C/Cu is of type D, its inverse is a tridiagonal Z- 
matrix with all row sums zero except the first, which corresponds to the (m + l)-st row 
of ^. 



Now one can easily get the entries of (C/C22) ^ and (C/Cu) ' using Theorem 4.5 
Furthermore, it follows that 



-biCu'Ei2{0/Cn)-' 



'±emeJ(C/Cu)-' 
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O-mO'tn+l — blb2 

Similarly, 

-b2C^2^E2iiC/C22)-^ = r^eie„_ 



^2 T 



This establishes the result in Case I. 



Case II {am = or Om+i =0): if am = or am+i = (or both), form a new cyclops from 
C by replacing the m-th or the (m + l)-st parameter (or both) by real numbers that 
approach zero. The result then follows from Case I and continuity. ^ 

Notice that if C is as in the previous theorem, then all row sums and column sums of ^ = 
are zero, except at least one of the m-th or the (m + l)-st (for otherwise C would be singular). 



Corollary 4.7 Let C € IR"^^'" be a nonsingular cyclops with eye m+ and parameters 0^,61,62. 
Then is a Z-matrix if and only if the following two conditions hold: 

(i) 61 = or sgn(6i) = sgn(amam+i - bib2), 

(ii) 52 = or sgn(62) = sgn(amam+i - bib2). 



Theorem 4.8 Let C G H^'" be a nonsingular cyclops with eye m+ and parameters ai,bi,b2, 
whose inverse is a Z-matrix. Let 

X = { k — j \ ajttk — bib2 < 0, j < m, k > m + 1 }, 

T = { k — j \ ajak — 6162 > 0, j < m, k > m + 1 }, 



and define 



X = 



min(x) X 7^ 
n + 1 otherwise, 



y 



min(T) ifx^^ 
n + 1 otherwise, 

r = number of positive aj with j < m, 
t = number of positive aj with j > m. 

Then the following hold: 

(i) If detC > and am > 0, then C E L^ (i.e., C is an inverse M-matrix.) 
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(ii) //detC > and < 0, then C £ L^, where 

s = n — 1 — max{m, n — m,n — x + 1}. 

(iii) //detC < 0, then C G Ls, where 

s = n — 1 — max{r, n — y + 1}. 



Proof: Let i? be any principal submatrix of C, partitioned as in Then i? is one of three 

types: 



1. i? is a principal submatrix of Cn, in which case, by (4.5), deti? has the same sign as the 
parameter Oj with the largest index contained in B. 

2. i? is a principal submatrix of C22, in which case det-B has the same sign as the parameter 
Oj with the smallest index contained in B. 

3. B is neither a principal submatrix of Cu nor of C22] in this case deti? has the same sign 
as ajak — ^1625 where aj has the largest index less than m contained in B, and has 
the smallest index greater than m + 1 contained in B. 



If detC > 0, by Corollary iA we have that amdm-i — ^1^2 > 0. Since C is an inverse Z-matrix, 
it follows from our previous results that 61 > 0, 62 > 0. Hence am and am+i are nonzero 
and have the same sign. If am > 0, it follows that C is nonnegative and hence an inverse 
M-matrix, i.e., (i) holds. If am < 0, then am+i < and detCn, detC22 are both negative; 
thus C has negative principal minors of sizes m x m and (n — m) x (n — m). We need also 
consider submatrices of the third type; the largest such submatrix with a negative determinant 
is of size (n — x + 1) x (n — x + 1). By Theorem 2.1 applied to A = C~^, we have that (ii) 
holds. 

If detC < 0, then s is determined by the size of the largest submatrix of C with a positive 
determinant. The largest submatrix of Cu with a positive determinant is r x r. The largest 
submatrix of C22 with a positive determinant is t x t. The largest submatrix of C of the third 
type is of size (n — y + 1) x (n — y + 1), and (iii) follows. ^ 



Example 4.9 The following example illustrates a cyclops and its tridiagonal inverse, com- 



puted by Theorem 4.6 



C 



4 


3 


2 


-1 


-1 


-1 


3 


3 


2 


-1 


-1 


-1 


2 


2 


2 


-1 


-1 


-1 


-4 


-4 


-4 


1 


1 


1 


-4 


-4 


-4 


1 


2 


2 


-4 


-4 


-4 


1 


2 


3 



1 -1 

-1 2 -1 
-1 0.5 
-2 



-0.5 

-1 
-1 2 -1 
-1 1 



14 



Note that, as shown in Theorem 4.6, all row sums and column sums of are zero expect 
the 3~rd and the 4-th. Moreover, using Theorem 4.3, one easily obtains that detC = —2. 
Applying Theorem 4.8, we have y = 3, r = 3, t = 3, and thus C^^ G Li. 



Next, we will characterize generalized ultrametric matrices whose inverses are tridiagonal. We 
begin with the irreducible case. We remind the reader that C G TR^'^ is called irreducible if its 
directed graph, r(C), is strongly connected (see e.g., Also recall that C is called totally 
nonnegative if all its minors are nonnegative. 



Theorem 4.10 Let C G TR^'^ be a nonsingular matrix. Then the following are equivalent: 

(i) There is m & (n) such that is an irreducible row and column diagonally dominant 

tridiagonal M-matrix whose row and columns sums are all zero, except at least one of 
the m-th or [m + \)-st. 

(ii) There is m ^ (n) and Oj, 6i, 62 G H such that C is a cyclops with eye m+ and parameters 

Oj , 61 , 62 satisfying 

min{am,am+i} > max{6i,&2} and min{6i,62} > 0. 

(iii) C is an irreducible generalized ultrametric matrix whose inverse is tridiagonal. 

(iv) C is a totally nonnegative irreducible generalized ultrametric matrix. 



Proof: Let C = [dj] and A = C-^. 

(i) implies (ii): By [^, Theorem 3.2] applied to A and , cu = Cij = cji for alll < j < i < m 
and Cii = Cjfc = c^j for all m + 1 < i < /c < n. Since A is a tridiagonal M-matrix, by |lO| , 
Theorem 4.1], c„fc = 'Ttrt/'' = ''rn,m+i and c,k = ^"f^ = c^k for ah j < m, and 
k > m + 1. Similarly, Cjk = Cm+i,k = Cm+i,m for all j > m + 1, fc < m. Thus C is a cyclops with 
eye m+ and parameters Oj = ca, bi = Cm,m+i, &2 = Cm+i,m- Since j4 is a nonsingular M-matrix, 
all principal minors of C are positive and thus amdm+i — fei&2 > 0. Since A is irreducible, the 
entries of C are all positive. Lastly, the inequality min{am, 0^+1} > max{bi,b2} follows from 
the facts that C is an entrywise positive cyclops and that is a row and column diagonally 
dominant matrix (and thus each diagonal entry of C is greater than or equal to the other 
entries in the corresponding row and column, see |^, Theorem (3,5)]). 

(ii) implies (iii): Notice first that the conditions on the parameters of the cyclops imply that 



C is an generalized ultrametric matrix. By Corollary 4.7, C is invertible, and by |10| , Theorem 
4.1], is tridiagonal. 

(iii) implies (i): We only need to show that the row and column sums are as claimed. 



Claim I: For all j < i < k, either ca = Cik = c^i or ca = Cij = Cji. 
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Proof of Claim I: ^From [l^, Theorem 4.1], we see that for all j < i < k, cjk = ^i^^ and 
Ckj = ^ e''^'^ • Consider the triangle on vertices i,j,k (see Definition 2.3]). If i is the 
preferred vertex, then cji = c^i and Cij = Cik, and hence cjk = c^j and Cjk > max{cji,Cik}- 
But since C is a generalized ultrametric matrix, cu > max{cji, Cj^}, thus cu = Cij = Cji = 
Cjk = Ckj = Cik = Cki- If j is preferred, it follows that cjk = cji and Ckj = Cij and thus 
Cu = Cik = Cki- If k is preferred then cu = Cij = Cji. This establishes Claim I. 

Claim II: If p is the first nonzero row sum of A, then all other row sums are zero, except 
possibly the {p + l)-st. 

Proof of Claim II: For all j < p + 1, Theorem 3.2] implies Cp+i^p+i ^ Cj^p+i. Suppose 
there is q > p + 1 such that q-th. row sum is nonzero. Then for all k > p + 1, Theorem 
3.2] implies Cp+i^p+i / Ck^p+i- If we apply Claim I with i = p + 1 we have a contradiction 
that establishes Claim II. 

Claim III: If the q-th row sum of A is nonzero, then the column sums of columns 1,2,..., q — 2 
are zero. 

Proof of Claim III: By ||9|, Theorem 3.2], for all i < q and k > q, cu 7^ Cki and hence by 
Claim I, for all j < i, cu = Cij = Cji. If we now apply Theorem 3.2] to , we see 
that columns 1, 2, . . . , g — 1 must have zero column sums. This establishes Claim III. 

Suppose now that the p-th row sum of A is the first nonzero row sum. Let P be the permutation 
matrix which reverses the order of the indices 1, 2, . . . , n. 

If the {p + l)-st row sum of A is also nonzero, then by applying Claim III to A and PAP^ , 
with q = p and with q = p + 1, we see that the only possible nonzero column sums are the 
p-th and the {p + l)-st. Taking m = p, the implication is proven. 

If the {p + l)-st row sum is zero, then by applying Claim III to A and PAP^, with q = p, we 
see that the only possible nonzero column sums are the {p — l)-st, the p-th and the (p+ l)-st. 
By Claim II applied to A^, either the (p — l)-st or the (p + l)-st sum is zero. By choosing m 
appropriately to be either p — 1 or p, the implication is proven. 

(iii) if and only if (iv): Follows from the results in Q or [13|, and in g 

Corollary 4.11 A matrix A is of type D with parameter ai > if and only if A^^ is a 
tridiagonal M-matrix with the only nonzero row and column sums being the first. 

We say that C is a G-cyclops if it is nonsingular and satisfies any of the equivalent conditions 



of Theorem 4.1C. We also refer to a matrix all of whose entries are equal as a fiat matrix. 



Theorem 4.12 Let C be a nonsingular matrix that is reducible but not completely reducible. 
Then the following are equivalent: 
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(i) Either C or C'^ is of the form 



B :-- 



( Bii Bi2 Bi3 
B22 B23 



V 



Blm \ 
B2m 



Bm—l,m—l 





B. 



m—l,m 

B„ 



where Bn is a G-cyclops whose last column is of constant value given by fu; Bmm is a 
G-cy clops whose first row is of constant value given by fmm! each Bgt, 1 < s < t < m, 
is a fiat matrix whose fixed value is given by f^t; each Bgg, 2 < s < m — 1, is either a 
positive number fss or an entrywise positive matrix of the form 



B., 



fss fl 
9ss fl 



with Qss < fss', o-nd for some 2 < q < m — 1, 

/ll > /22 > . . . > fgq > 0, 
< /cj+l,cj+l < fg+2,q+2 < • • • < fmm, 

fit = f2t = ■■■ = ft-i,t = ftt for 1 <t <q, 
fss = fs,s+i = . . . = fs,m-i = fsm for q < s < m, 
< fst = fq,q+i < va.m.{fqq, for I < s < q < t < m. 

(ii) C is a generalized ultrametric matrix whose inverse is tridiagonal. 

(iii) C is a totally nonnegative generalized ultrametric matrix. 



Proof: 

(i) implies (ii): It is easy to see that the conditions on the fij guarantee that C is a generalized 
ultrametric matrix and that it satisfies []lO| , Theorem 4.1 (ii)]. 

(ii) implies (i): Let C = [cij] and A = C^^ = [oij]. Since A is a tridiagonal M-matrix, it 
must satisfy |[l^, Theorem 4.1 (ii)]. 

We begin by showing that C or must be block upper triangular with no zero entries in 
or above the diagonal blocks. Suppose that A has a zero entry on the superdiagonal and a 
zero entry on the subdiagonal. For simplicity, we will assume that ajj+i = and k > j is 
the smallest integer such that a^^i^k = (otherwise take A = {C'^)~^). If k = j then A is 
completely reducible contradicting our hypothesis, hence we will assume that j < k. Since 
j does not access j + 1 in T{C), by [15, Lemma 2.2], Cjj+i = 0. Similarly c^+i^k = 0. By 



Theorem 4.1 (ii)], cji = '^^+^'^^+^-' = for all i > j + 1, and Ck+i,i = Hfe±M£« = q for 



all / < k. Consider the triangle (see |1C, Definition 2.3]) on j,k,k + 1. Either Ckj = or 



Cfc,fc+i = 0. If Cfc = then by [15, Lemma 2.2] ak^k+i = and hence A is completely 
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reducible, contradicting our hypothesis. If c^j = 0, then by Lemma 2.2], k does not 
have access to j in T{A) and hence there must be an i with j < i < k such that aj+i,i = 0, 
contradicting the minimahty of k. So either the superdiagonal of or the subdiagonal of 
contains only nonzero entries. It follows that either C or C"^ is block upper triangular, 
as represented hy B. It remains to show that the blocks of B = [bij] are as claimed. Without 
loss of generality, assume that all the entries on the superdiagonal of A are nonzero. 

Since both B and A are block upper triangular, Bgs = {Ass)~^- Since Agg is irreducible, Bgg 



must be a G-cyclops by Theorem 4.10. Since the superdiagonal entries of A are nonzero, by 



the results in [^] each Bgt is an entry wise positive matrix for all s < t. 

Notice that if i < j and i and j are in different blocks of B then bji = 0. This fact will be used 
without further remark whenever triangles are considered in the remainder of this proof. We 
will also write j G s to mean that bjj is in the block Bgg- 

Let j G 1, G 1 and / ^ 1 with j < k. By considering the triangle on j,k,l, we see that 
bji = bki- But bji = "^^) hence bjk = bj-k and Bn is as claimed. A similar argument shows 
that Bmm is as claimed. 

Consider B„-, Bgs, and Btt with r < s < t. Let i G r, j ^ s, k & s, I & t. Consider 
the triangle on Then bu = min{bij,bji}. By Theorem 4.1], bu = and hence 

bjj = max{bij,bji}. Similarly, bkk = nT-0'X{bik,bki}. From the triangle on i,j,k we see that 
bij = bik < max{bjk, bkj}- From the triangle on j, k, I we have that bji = bki < max{bjk, bkj}- 
But then bjj = max{bij,bji} < max{bjk,bkj} < bjj. Hence equality must hold throughout. 
Using the corresponding inequalities for k, we can conclude that bjj = max{bjk,bkj} = bkk- If 
j < k, then bik = implies that fss = bjk = bjj > bkj = gss- Since Bss = {Ass)~^, B^g 

must be a nonsingular G-cyclops and hence can only be as claimed for 2 < s < m — 1. 

Let now h = min{j \ bjj = bjk for all k > j} {h is well defined since b^m — b^k fo^ ^ ^ j) 
Consider the r~th diagonal block so that h ^ r. By the equalities in the above paragraph 
if /i — 1 G r then = b^^i^h = b^h = b^k = ^/i-i.fc, contradicting the minimality of h 

Hence h — 1 ^ r. Set q = r — 1. If r < m, then by the choice of h and the triangle on h,h + l,k 
hhh = bh,h+i = bhk = min{bh,h+i, bh+i,k} for all > /i + 1, which implies that bh,h+i < hh+i,k 
Thus = 'max{bhhj^i,bhj^i^k} = ^h+i,fc for all /c > h + 1. We can now repeat this 

argument for /i + 2, /i + 3, . . . , up to largest index in the (m — l)~st diagonal block to conclude 
that for all g < s < t with j G s and k t, fss = bjj = bjk = fst and thus the Bst are as claimed. 
For any j < h, with j ^ 1, by the choice of h and the inequalities in the above paragraph, 
we see that bjj = bij for all i < j. Hence Bst must be as claimed for all s < t < q. For 
s < g < t, let j G s and G t. Then bjk = = bjh and bjh = ^^trfcr^- considering 

the triangle on j,h — l,/i, we have that bjh = min{bj^h-iibh-i,h} = fnin{bh-i,h-i-,bh-i,h\ = 
bh^i,h < min{bh-i,h-i, bhh}- Thus fst = bjk = bh^i,h = fq,q+i- 

(ii) if and only if (iii): Follows from the results in B or [13|, and in Q. g 



Example 4.13 The following matrix illustrates a matrix that satisfies the conditions of The- 
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orem 4.12| (and hence it is a reducible, totally nonnegative, generalized ultrametric matrix 
whose inverse is tridiagonal). 



C 



1 9 


1 1 

L X 


1 n 


Q 


Q 


7 


o 


(J 


(J 


o 


o 


o 


11 


1 -1 

11 


10 


9 


9 


'-T 

7 


5 


5 


5 


5 


5 


5 


10 


10 


10 


9 


9 


7 


5 


5 


5 


5 


5 


5 











9 


9 


7 


5 


5 


5 


5 


5 


5 











8 


9 


7 


5 


5 


5 


5 


5 


5 

















7 


5 


5 


5 


5 


5 


5 




















6 


6 


6 


6 


6 


6 























7 


7 


7 


7 


7 























6 


7 


7 


7 


7 





























8 


8 


8 





























8 


9 


9 





























8 


9 


10 



Finally, Theorems 4.10| and 4.12 yield the following result. 



Theorem 4.14 Let C G TR^'^ be a nonsingular matrix. Then the following are equivalent: 



(i) C is the direct sum of matrices of the forms given in Theorem 4.1(\ (ii) and Theorem 4-l^ (i)- 

(ii) C is a generalized ultrametric matrix whose inverse is tridiagonal. 

(iii) C is a totally nonnegative generalized ultrametric matrix. 
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